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Recent numerical simulations of dynamo action resulting from rotating convection have revealed 
some serious problems in applying the standard picture of mean field electrodynamics at high 
values of the magnetic Reynolds number, and have thereby underlined the difficulties in large- 
scale magnetic field generation in this regime. Here we consider kinematic dynamo processes 
in a rotating convective layer of Boussinesq fluid with the additional influence of a large-scale 
horizontal velocity shear. Incorporating the shear flow enhances the dynamo growth rate and also 
leads to the generation of significant magnetic fields on large scales. By the technique of spectral 
filtering, we analyse the modes in the velocity that are principally responsible for dynamo action, 
and show that the magnetic field resulting from the full flow relies crucially on a range of scales 
in the velocity field. Filtering the flow to provide a true separation of scales between the shear 
and the convective flow also leads to dynamo action; however, the magnetic field in this case has 
a very different structure from that generated by the full velocity field. We also show that the 
nature of the dynamo action is broadly similar irrespective of whether the flow in the absence of 
shear can support dynamo action. 



Oh 



Key words: 



1. Introduction 

One of the outstanding theoretical problems in astrophysical MHD is to account for the gener- 
ation of global scale magnetic fields, as detected in many cosmic bodies. These are generally held 
to be produced by a hydromagnetic dynamo process, in which the magnetic fields are maintained 
against resistive effects by induction due to the plasma motions. It is, however, hard to provide a 
convincing theoretical explanation of how such large-scale fields — i.e. those with a significant 
component on scales much larger than that of the plasma motions responsible for their generation 
— can be generated. 

The traditional theoretical approach to explaining th e generation of large-scale magnetic field s 
is via mean field electrodynamics (see, for example, lMoffattl[l978l: iKrause & Radler1ll98Qh . 



Here, the evolution of the mean (large-scale) field is governed by the mean induction equation, 

-^ = Vx([/ xB ) + Vx5 + r/V 2 S , (LI) 

where Bq represents the mean magnetic field, Uo the mean velocity, £ the mean electromotive 
force (emf) and r] the magnetic diffusivity. The mean emf is defined by 

S = (u x 6), (1.2) 

where u and b represent the (small-scale) fluctuating velocity and magnetic fields, and angle 
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brackets denote a spatial average over intermediate scales. The theory proceeds on the assump- 
tion that £ is a linear functional of So, which leads to the expansion 

9B 0j j 



Si 



ctijB, 



0] 



■A 



ijk 



dx k 



(1.3) 



(Hug hes & Proctor I d2010h discuss the implications of a more general expansion procedure in- 
volving also temporal derivatives of the mean field.) In the kinematic regime, in which the field 
is assumed to exert no back-reaction on the flow, the components and fiijk depend solely on 
the properties of the velocity field and on the magnetic diffusivity. The symmetric part of the a 
tensor (the so-called 'a-effect') leads to field amplification, and can be non-zero only in flows 
that lack reflectional symmetry, such as helical flows. For isotropic turbulence, = /3e^ 
and the scalar /3 can then be identified as a turbulent diffusivity; in general though, fajk has a 
much more complicated interpretation (see iKrause & Radlenll980l) . In most astrophysical ap- 
plications, the mean field is considered to be axisymmetric; it can then be written, in cylindrical 
polar coordinates (s, 0, z), as Bq = V x (A(s, z)e < f ) ) + B(s, z)e ( / ) . Under the strongest sim- 
plifying assumptions of isotropic turbulence and an azimuthal mean flow, the mean field is then 
described by the following two equations: 

?A=aB + fj(s,z)(v 2 -^)A, (1.4) 

-=s(B P -V) W + (Vx (aB P )) • e + fj f V 2 - ^ J B + - Vfj • V(sB), (1.5) 

where Bp denotes the poloidal field, scj(s, z)e ( j ) is the mean flow and fj(s, z) = 77 + f3 is the 
total magnetic diffusivity. It is necessary for dynamo action that the coupling terms in these 
equations are non-zero; the dynamo cycle depends crucially on the generation of poloidal field 
from toroidal, and, conversely, the generation of toroidal field from poloidal. The former requires 
the a-effect, whilst the latter can result from either the a-effect or from shearing of poloidal field 
by the differential rotation, the cj-effect. The resulting dynamos are designated, respectively, as 
a 2 or auj -dynamos. Closed form expressions for a and f3 can be obtained only under simplifying 
assumptions, notably small values of the magnetic Reynolds number on the fluctuating scale, 
or short correlation times for the flow. Neither of these applies, however, in the astrophysical 
context. In consequence, astrophysical modelling typically involves adopting plausible, albeit 
arbitrary, spatial forms and amplitudes of a and f3. 

Recent research has focused on attempts to measure the a-effect directly in numerical sim- 
ulations of turbulent flows, either in forced helical tur bulence (e.g. Cattaneo & Hughes I ll996t 



Brandenburg 2001; Cattaneo, Hushes & Thelen 2002) or in rotating turbulent convection (e.g. 
Cat taneo & H ughes 2006; Hug hes & Cattaneo! 2008k . Mean field coefficients can only be prop- 



erly determined if there is adequate separation between the small scales of the turbulence and 
the system size; for convective turbulence, this is most readily accomplished by adopting the 
relatively simple system of plane layer, Boussinesq convection. The study of large-scale dy- 
namo action in this s ystem has quite a l ong h istory, dating back to the pioneering papers of 
IChildress & Sowardl ( 1972 ) and lSowardl d 19741). Subsequently, there have been a number of nu- 



merical investigations of the probl em (e.g. lSt. Pierre 119931: Jones & Roberts hood iRotvig & Jones 



12002k Istellmach & Hansen 1120041) . These have considered magnetic field generation in domains 
with an 0(1) aspect ratio, driven by mildly supercritical convection at fairly rapid rotation rates; 
the resulting dynamo can then be interpreted as a mean field a 2 dynamo. 

More recently, with ever-improving computational perf ormance, it has become possible to in- 

vestigate more turbul ent regimes at larger a s pect ra tios fe.g. lCattaneo & Hughes 120061 iKapyla. Korpi & Brande nburg 
2010|) . The paper of ICattaneo & Hughes I (I2006|) suggested two significant problems with the 



standard mean field picture. The first is that when the convection is sufficiently vigorous, it acts 
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as a small-scale dynamo, despite the flow being significantly helical. There is no tendency to 
generate large-scale field; indeed, the spectrum of magnetic energy is essentially ident ical to that 



resulting from the small-scale dy namo generated by turbulent non-rotating convection (ICattaneo 



19991: ICattaneo & Hughes 1 l2006h . Attempts to measure the a-effect by imposing a uniform hor- 



izontal field for turbulent convection just below the dynamo threshold, but still at high magnetic 
Reynolds number Rm, reveal the second problem. Despite averages being taken over many con- 
vective cells, the a-effect exhibits significant temporal variations about a mean value that is much 
smaller than the characteristic speeds of the flow. Surprisingly, a coherent helicity distribution 
does not lead to a significant a-effect. 

However, rather different conclusions were reported for the convective dynamo simulations 
of iKapyla. Korpi & Bra ndenburg 1 d2010h . who obtained significant mean fields and sizeable a- 
effects. The differences can be attributed to a number of factors: (a) the a-effect depends on 
the horizontal correlations of the turbulence — which, in turn, depend on the degree of super- 
criticality and the rotation rate; (b) different calculations employ different magnetic boundary 
conditions; (c) there are different definiti ons of the a-e ffect and the means of measuring it — 
e.g. th e traditional imposed field method (Mo ffattlll978h. the test field method ( Schrinner et al. I 



20071) and the 'resettin g' method (jOs sendriiver et al. 2002). A detailed discussion of all of these 
issues can be found in lHughes. Proctor & Cattaneo ( 201 lb 



Although the precise nature of the dynamo mechanism in these simulations remains uncertain, 
there is no doubt that large-scale fields are observed in nature and it is therefore important to 
identify other mechanisms that may lead to such fields. Since most astrophysical bodies possess 
a large-scale differential rotation, it is natural to incorporate a large-scale shear flow into the con- 
vection model and explore the consequences for any dynamo act i on. In this paper we carry out 



such a programme, building on the results of iHughes & Proctor I (120091) and lProctor & H ughes 



(1201 II) . who were the first to show that a combination of small-scale convection and large-scale 
velocity shear could lead to magnetic field growth on large scales. In order to obtain scale sepa- 
ration between the convection and the shear flow, we consider a uni-directional horizontal flow, 
dependent only on the other horizontal direction. We concentrate in this paper solely on the kine- 
matic dynamo problem, in which the back-reaction of the magnetic field on the velocity via the 
Lorentz force is ignored; thus we examine in some detail the nature of the generation mechanism, 
but we do not address the means by which magnetic field growth is saturated. 

The paper is organized as follows. In §[2] we discuss the various ways in which a large-scale 
shear flow may affect the dynamo process. The mathematical formulation of the problem is 
contained in §[3] In §[4] we consider the introduction of a velocity shear into a convective flow 
that, in the absence of shear, does not act as a dynamo; we describe first the characteristics of the 
flow and then those of the magnetic field that it generates. Section [5] looks in more detail at the 
dynamo process, through considering 'filtered' flows, in which only certain scales in the velocity 
are retained. In §[6] we consider a more vigorous convective state than in §ffl one that supports 
dynamo action even in the absence of shear, in order to determine whether this is a crucial factor 
in the nature of the ensuing dynamo action. In the concluding §|7]we discuss the implications of 
our results and their relation to parallel studies of dynamos driven by a combination of forced 
turbulence and uniform shear. 



2. The Influence of Velocity Shear on Convective Dynamos 

Before describing our results it is instructive to consider the various possible ways in which a 
large-scale velocity shear may influence the nature of dynamo action driven by rotating convec- 
tion. A number of possibilities suggest themselves. 

At high Rm, rotating convective turbulence, in the absence of shear, can induce large local 
emfs, but these are decorrelated in space and time, leading to a small net a-effect. In conse- 
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quenc e, any dynamo field generated in extended domains is predominantly small-scale (ICattaneo & Hughes 



2006h . It is though conceivable that a coherent large-scale shear may impose more order on the 



correlations and, in so doing, enhance the a-effect. Alternatively, even if the mean emfs remain 
very small, a large shear may be able to compensate for a feeble a-effect (or a more complicated 
mean-field process) to make a viable two-scale dynamo; in the classical mean field picture it is 
the product of a and uj that controls the efficiency of the dynamo. 

A rather different possibility is that enhanced dynamo action may depend on the interaction of 
a wide range of scales, from the largest scale of the shear to the convective cell size. In an extreme 
version of this, dynamo action might result solely from the interactions between the large-scale 
shear and the induced motions on a similarly large scale; this would then be effectively a small- 
scale (i.e. one- scale) dynamo, but on the scale of the shear flow rather than that of the convection. 
We shall interpret our findings with these possibilities in mind, considering cases for which the 
convection in the absence of shear does, and does not, act as a dynamo. 

3. Formulation 

Following Cattaneo & Hughes l d2006h and Hughes & Cattaneo ( 2008), we consider thermally 



driven convection in a three-dimensional, Cartesian layer (0 < x, y < \d, < z < d) of Boussi- 
nesq fluid rotating about the vertical. The layer has angular velocity 17, density p, kinematic 
viscosity v, thermal diffusivity k and magnetic diffusivity 77. This basic model is then extended 
by the inclusion of a horizontal flow of the (dimensional) form 

U = U f(y/d)x , where f(y/d) = cos — , (3.1) 

where the total velocity is now u + Uq ; this is accomplished by replacing u with u + Uq in 
the governing equations except for the viscous term (equivalent to forcing the flow via the mo- 
mentum equation, but eliminating viscous transients). We adopt a periodic flow for con sistency 
with the periodic horizontal boundary conditions adopted in lCattaneo & Hughes I (120061) . For the 
purposes of this paper we shall restrict attention to kinematic dynamo action, so that the back- 
reaction of the Lorentz forces on the convection is neglected, as is appropriate for very weak 
fields. 

Following standard practice, we adopt the layer depth d, the thermal relaxation time d 2 /K, 
and the temperature drop across the layer AT as the units of length, time, and temperature re- 
spectively. All velocities are scaled with njd\ in particular, Uq below is now dimensionless. The 
governing non-dimensional equations for the velocity u, temperature perturbation and mag- 
netic field B can then be expressed as 

(d t -aV 2 )u + u-\7u + U (f(y)d x u + f f (y)uyX ) + aTa 1/2 z x u = -Vp + crRaOz , (3.2) 

(d t - (V 2 )B + u • VB + U f(y)d x B = BVu + U f'(y)B y x , (3.3) 

(d t - V 2 )6 + u • W + U o f(y)d x = u z, (3.4) 

V • B = V • u = 0, (3.5) 
where w is the vertical velocity, and denote s the temperature fluctuations relative to a linear 



background profile (e.g. IChandrasekar 1 1 1 9611) . As noted in the introduction, we consider here 
only the kinematic dynamo problem, and thus the Lorentz force is omitted in the momentum 
equation (13.2k the problem is then linear in the magnetic field, the scaling of which is arbi- 
trary. Five dimensionless parameters appear explicitly in the governing equations: the Rayleigh 
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number Ra = gaf3d A / ' nv (where g is the gravitational acceleration, a is the coefficient of ther- 
mal expansion and /3 is the superadiabatic temperature gradient), which measures the strength 
of thermal buoyancy relative to dissipation; the Taylor number Ta = AVt 2 d 4 / v 2 ; the kinetic and 
magnetic Prandtl numbers 

v v 
a = — and a rn = — , (3.6) 

K T] 

and the dimensionless speed Uq. Additionally there is the choice of the aspect ratio A. 

The purely hydrodynamic solution is evolved until it reaches a stationary state, starting from 
an initial condition of a small perturbation to the shear flow (13.11) . We then consider the dynamo 
action resulting from such stationary states. It should be noted that although a flow with a large- 
scale component (i.e. with the same spatial dependence as the 'target flow' (13.1b ) does indeed 
occur, its amplitude may differ appreciably from Uq\ the hydrodynamic state that ensues depends 
on interactions between the shear flow and convection and, possibly, on instabilities of the shear 
flow itself. We also introduce the derived quantity 

S = 27rU £/\u rms , (3.7) 

where £ and u rms are estimates, respectively, of the horizontal scale of the convection and of the 
typical velocity in the absence of shear; S provides a measure of the competition between shear 
and convection. 

In the horizontal directions we assume that all fields are periodic with periodicity A. In the 
vertical we consider standard illustrative boundary conditions on the temperature and velocity 
fields, namely that the boundaries are perfect thermal conductors, impermeable and stress-free. 
Formally these correspond to 

= w = d z u = d z v = at z = 0, 1. (3.8) 

The natural average in this system is one over horizontal planes, which involves averaging over 
many convective cells. From the point of view of generating large-scale fields with the simplest 
vertical structure, it is therefore preferable to choose perfectly conducting boundary conditions, 
for which the field is purely horizontal, thereby admitting field configurations with only one node 
in the vertical. Thus we choose 

B z = d z B x = d z B y = at z = 0, 1. (3.9) 

Equations (13.21) - (13.51) are solved numerically by standard pseudo-spectral methods optimized 
fo r machines with parallel architectur e. Details concerning the numerical methods can be found 
inlCattaneo. Emonet & Weiss I d2003l). 



ICattaneo & Hughe7 l d2006h and lHughes & Cattaneo I (1200 8) explored dynamo action and mean 



emf generation in systems with the fixed values of Ta = 500 000, a = 1, a m = 5 and for values 
of Ra between 80 000 (slig htly above the onset of conv ectio n) and 1 000 000; the onset of dynamo 
action is at Ra w 170 OOO. lHughes & Proctorl d2009h and lProctor & Hughes I d201lh examined 



the influence of a range of values of shear amplitude Uo for the same fixed values of Ta, a and 
cr m , for Ra = 150 000 (for which there is no dynamo in the absence of shear) and aspect ratio 
A = 5. Here we consider extended spatial domains, with A = 10 and, for a few runs, A = 20, 
and consider both Ra = 150 000 and Ra = 250 000 (for which there is a small-scale dynamo in 
the absence of shear). For Ra = 150 000, ^ rms w 60 (in the absence of shear), and the width of 
the convective cells is comparable with, though a little smaller, than the layer depth; taking I w d 
leads to S ~ Uo/600 in this case. In 2] and ^5]we concentrate on the case of Ra = 150 000; 
similarities and differences for the case of Ra = 250 000 are discussed in The numerical 
resolution and parameter values for all the simulations presented in this paper are summarized in 
Table[Q 
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Ra 




Ta 


Uo 


a 




A 


N x xN y x N z 


1 


5 


X 


10 5 


5 


X 


10 5 





1 


5 


10 


512 


X 


512 x 97 


1 


5 


X 


10 5 


5 


X 


10 5 


100 


1 


5 


10 


512 


X 


512 x 97 


1 


5 


X 


10 5 


5 


X 


10 5 


200 


1 


5 


10 


512 


X 


512 x 97 


1 


5 


X 


10 5 


5 


X 


10 5 


300 


1 


5 


10 


512 


X 


512 x 97 


1 


5 


X 


10 5 


5 


X 


10 5 


400 


1 


5 


10 


512 


X 


512 x 97 


1 


5 


X 


10 5 


5 


X 


10 5 


500 


1 


5 


10 


512 


X 


512 x 97 


1 


5 


X 


10 5 


5 


X 


10 5 


600 


1 


5 


10 


512 


X 


512 x 97 


1 


5 


X 


10 5 


5 


X 


10 5 


700 


1 


5 


10 


512 


X 


512 x 97 


1 


5 


X 


10 5 


5 


X 


10 5 


800 


1 


5 


10 


512 


X 


512 x 97 


1 


5 


X 


10 5 


5 


X 


10 5 


900 


1 


5 


10 


512 


X 


512 x 97 


1 


5 


X 


10 5 


5 


X 


10 5 


1000 


1 


5 


10 


512 


X 


512 x 97 


1 


5 


X 


10 5 


5 


X 


10 5 


2000 


1 


5 


10 


512 


X 


512 x 97 


1 


5 


X 


10 5 


5 


X 


10 5 


1000 


1 


5 


20 


1024 


X 


1024 x 97 


1 


5 


X 


10 5 


5 


X 


10 5 


2000 


1 


5 


20 


1024 


X 


1024 x 97 


2 


5 


X 


10 5 


5 


X 


10 5 





1 


5 


10 


512 


X 


512 x 97 


2 


5 


X 


10 5 


5 


X 


10 5 


200 


1 


5 


10 


512 


X 


512 x 97 


2 


5 


X 


10 5 


5 


X 


10 5 


400 


1 


5 


10 


512 


X 


512 x 97 


2 


5 


X 


10 5 


5 


X 


10 5 


600 


1 


5 


10 


512 


X 


512 x 97 


2 


5 


X 


10 5 


5 


X 


10 5 


800 


1 


5 


10 


512 


X 


512 x 97 


2 


5 


X 


10 5 


5 


X 


10 5 


1000 


1 


5 


10 


512 


X 


512 x 97 


2 


5 


X 


10 5 


5 


X 


10 5 


1200 


1 


5 


10 


512 


X 


512 x 97 


2 


5 


X 


10 5 


5 


X 


10 5 


1400 


1 


5 


10 


512 


X 


512 x 97 


2 


5 


X 


10 5 


5 


X 


10 5 


1630 


1 


5 


10 


512 


X 


512 x 97 



Table 1 . Summary of the parameter values and numerical resolution for the simulations. 



4. Flows and Fields in Rotating Sheared Convection 

4.1. Influence of Shear on Convection 

Understanding the interactions between a shear flow and rotating convection is a complex prob- 
lem, of relevance for both stellar and planetary physics. The majority of work has focused on the 
case of a shear flow dependent on the vertical direction, arising from relat ive motion of the hori- 



zontal boundaries or, alternatively, from a fictitious force (see, for examp le, Hathawav, Toomr e & Gilman 
1980: lHathawav & Somervillell983[[T986HKropp & Bus sell 99 ll: [Matthews & Cox 111 9971: ICox*~ 



1998b . The hydrodynamical problem of a horizontally dependent shear flow, the case we consider 



here, may also be of relevance in planetary atmospheres, and has been examined in the nonlinear 



regime by lHathaway & Somervillel (1987). 



On increasing the amplitude Uo of the target shear flow, various regimes can be identified 
in the nature of the resulting convection. These are demonstrated in figure [TJ which plots the 
ratio of the kinetic energy in the target flow to the total kinetic energy as a function of time for 
three different values of Uo, and figure [2 which shows the corresponding density plots of the 
temperature fluctuations close to the upper boundary, together with that of the non-sheared state. 
For Uo = 300 (shear parameter S ~ 0.5) the convection is such as to decrease the energy in the 
shear mode from its target value; note from figure [Ha) that, at least for this value of Uo, a long 
temporal integration is needed in order to determine the final stationary state. In figure 06) it 
is possible to detect a large-scale vortex underlying the small-scale convection. For 0(1) values 
of S, the kinetic energy in the shear flow is comparable with its target value, and this mode 
dominates the total kinetic energy (e.g. figure \Hb)) (It should be noted that the hydrodynamic 
state has been evolved for much longer than shown in figureQJfr), wim no transition to a different 
state.) As shown in figure Etc), the shear leads to a clear elongation of the convective cells, 
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f , , , i , , i , , , i , , , i , , , 1 9R0 L , , , , i , , , , i , , , , i , , , , 1 OO f,,,,,, , , , i , , , , i , , , , 1 
2 4 6 8 10 0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 

t t t 



FIGURE 1 . Temporal evolution of the energy in the target shear flow mode (i.e. the mode with wavenumbers 
k y = 1, k x = k z = 0) normalized by the total kinetic energy of the flow: (a) Uo = 300, (b) Uo = 1000, 
(c) Uo = 2000. 




FIGURE 2. Snapshots of the temperature perturbations close to the upper boundary for Ra — 150 000 and 
four different values of the shear flow: (a) U = 0, (b) U = 300, (c) U = 1000, (d) U = 2000. 
White denotes hot fluid, black cool fluid. The target shear flow is Uo — Uo cos(27ry/A)x , where x is the 
direction of the abscissa. 
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together with significant inhomogeneity between the two halves of the domain in the y direction. 
For < y < A/ 2, the vorticity augments the underlying vorticity due to the rotation of the layer, 
whereas for A/2 < y < A it tends to reduce it. The net underlying vorticity in the z-direction can 
be expressed in dimensionless form as 

1/2 , 2tt/7 . 2th/ 

Ta 1 H — sin—. (4.1) 

A A 

Clearly (when Uo is positive) the underlying vorticity has the smallest absolute value when y = 
3 A/4. The vorticity dynamics in the two halves of the layer is similar if Uo is very small or 
large; the maximum disparity between the two halves of the layer (in y) occurs when Uo ~ 
XTa 1 ! 2 /2tt & 1125 here. In figure[2tc), Uo is close to this optimal value, and it can be seen that 
convection is indeed most vigorous in the neighbourhood of y = 3 A/4. For Uo < 500 (S < 1), 
convection dominates in the sense that there are no streamlines extending across the domain. For 
larger values of Uo (e.g. figure[2tc)) a clear 'channel flow' is established in < y < A/2. For a 
range of values of the shear amplitude Uo, this shear-dominated flow remains stable. However, 
at yet larger values of Uo, the shear becomes unstable and the resulting flow reverts to being less 
shear-dominated, as can be seen by figure [He) for Uo = 2000. At these larger values of Uo, a 
large coherent vortex forms and the flow has a very different structure, with the convective cells 
expelled from the vortex (see figure Wid)). In this paper we shall concentrate principally on the 
nature of the dynamo action resulting from values of Uo for which the convection and shear flow 
can co-exist (e.g. figure Efc)) and for which there is a clear separation in their spatial scales. 
Figure [3] gives a measure of the planarity V for the flows with Uo = and Uo = 1000, where V 
is defined as the ratio of the horizontal to total kinetic energies, 

(U 2 ) 

n V ,z) = ^, (4.2) 

with angle brackets denoting an average over x. It can be seen that for Uo — 1000 the flow is es- 
sentially two-dimensional for much of the domain, with patches of fully three-dimensional flows 
centred around the turning points in the target shear flow. For Boussi nesq convection, the helicity 



distribution is anti- symme tric about the mid-plane (see, for example Jchildress & Soward 

E2Z2 



Cattaneo & Hughes 1120061) . However, the introduction of a shear flow in a rotating frame allows 
for differences between the domains y < A/ 2 and y > A/ 2; this is illustrated by figure [U which 
shows the relative helicity h(z) for the two halves of the ^-domain, where 

h(z) = ("-V*") (4 3) 

with the averages taken over horizontal planes. The helicity is significantly greater for y > A/ 2, 
where the background and flow vorticities are of the same sign. 

4.2. Kinematic Dynamo Action 

Figure [5] plots the dynamo growth rate as a function of Uo, for Ra = 150 000. It can be seen that 
the incorporation of velocity shear facilitates dynamo action, with the critical value of the shear 
amplitude given by Uo ~ 200 (i.e. S « 0.3). Further increases in Uo serve to enhance the growth 
rate, although there is no simple power law relationship. The levelling off in the growth rate for 
500 < Uq < 700 corresponds to a change in the nature of the flow regime, as described above. 
We have also calculated the dynamo growth rates for a domain that is twice as wide (A = 20); 
a comparison between the domains of differing sizes then requires replacing Uo with 2Uo, to 
keep the same value of S, according to definition (13. 71) . For the smaller value of S there is little 
influence of the box size, whereas at the larger value of S, although the growth rates are similar, 
there is a clear influence of the domain size, with the dynamo in the larger domain being more 
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Figure 4. Snapshots of h(z) (the horizontally averaged relative flow helicity) for y < A/2 (dashed line) 
and y > A/2 (solid line) for Uo = 1000. Exact antisymmetry about the midplane (z = 0.5) is recovered 
by time averaging. 



efficient. In any case, we do not necessarily expect close agreement, since although the shear S is 
the same for the A = 10, Uo = 1000 and A = 20, Uo = 2000 runs, the convective structures are 
independent of the box size whereas the region of, sa y, positive shear scales with A. The dynamo 
growth rates are also consistent with those found in Hugh es & Proctor in a domain of 

half the width (so the values of S in that paper should here be multiplied by two for comparison); 
the more extensive data that we now have makes it clear that the relationship between growth rate 
and shear is more complicated than the linear one with which the earlier data were consistent. 

Figure [6l which plots B x at the top of the layer, illustrates how the magnetic field changes 
with increasing Uo- At Uo = 200, essentially the smallest value of the shear flow that allows for 
dynamo action, although there is already some evidence of asymmetry between the two halves 
of the ^/-domain, there is still significant magnetic energy inO<?/<A/2(a 28% to 72% split in 
the energy of the B x field between the two halves of the y domain). At Uq = 300, the underlying 
vortex depicted in figure [26) has a clear influence, introducing a strong large-scale variation in 
magnetic field in the x-direction. With a further increase in Uo to Uo = 400, the flow enters 
the regime of 0(1) values of S, and the field shows a clear asymmetry between the two halves 
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Figure 5. Dynamo growth rates versus Uo forRa = 150 000 and A = 10 (asterisks). The diamonds show 
the growth rates when the shear amplitude is 2Uo with A = 20. 




FIGURE 6. Density plots of B x at the upper boundary, for increasing values of Uo'. (a) Uo — 200, (b) 
Uo = 300, (c) Uo = 400. Each plot is scaled individually; colour table as in figure [7] below. 



of the ^-domain. The field persists in this form for a range of Uo, until the shear amplitude is 
sufficiently great as to trigger an instability (figure Wid)). 

Figure [7] shows snapshots of B x and B y in the x?/-plane at the upper boundary, for Uo = 500 
and Uo = 1000. A movie of such plots (included as supplementary material) reveals clearly the 
advection of the magnetic field pattern by the velocity shear. Two important features can be noted. 
One is that the dynamo action is strongly inhomogeneous, being concentrated in A/2 < y < A. 
The other is that the stronger shear leads to pronounced stretching of the field structures in the 
x-direction. 

In order to obtain a quantitative description of the scales on which the field is being generated, 
two-dimensional Fourier transforms of B 2 and By are constructed, after both depth averaging 
and time averaging and having removed the exponential growth of the field. Figured] shows these 
plots for the range of horizontal wavenumbers ^ k x , k y < 32. The plots are scaled individually, 
but it should be pointed out that the bulk of the magnetic energy resides in B x , as is maybe to 
be expected from a flow that is strongly influenced by a shear U(y)x ; for the parameter values 
of figure^ {B 2 X )/{B 2 ) w 0.77, {B 2 y ) / {B 2 ) w 0.14, (B 2 )/{B 2 ) w 0.09 (where angle brackets 
denote an average over the fluid volume and time). The distribution over wavenumbers displays a 
marked asymmetry in k x and k y for (B 2 ) , but is roughly symmetric for (B 2 ) . For the former case, 
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Figure 7. Snapshots of the two components of the horizontal magnetic field at the upper boundary: (a) B x 
for U = 500, (b) B y for U = 500, (c) B x for U = 1000, (d) B y for U = 1000. The plots for each U 
are scaled individually; the colour bar ranges from — maxl^ | to max\B x \. 




FIGURE 8. Density plots of the time-averaged and depth- averaged Fourier transform of the magnetic energy 
in (a) B x and (b) B y for Uo = 1000. The horizontal and vertical axes are the x and y wavenumbers, 
respectively, in the range ^ k < 32. The plots are scaled individually, from blue (highest energy) to 
yellow (lowest). 
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Figure 9. One-dimensional spectra of the energy in (a) B x and (b) B y . The solid (dashed) lines show the 
spectra as a function of k x (k y ) after averaging over y (x). 



which provides the principal contribution to the overall magnetic energy, the dominant modes are 
k x = 1, k y = 1, 2, 3; the fall off with energy with increasing k x is significantly greater than that 
with increasing k y . An alternative representation of the distribution of the magnetic energy over 
wavenumbers is provided by figure |2 which shows the one-dimensional spectra of and By 
having summed over either k x or These spectra should be compared with that for the case 
of no velocity shear (e.g. figure 6 of Catt aneo & Hughes 1 l2006lt . in which the magnetic energy is 
peaked at the scale of the convective cells and falls off rapidly to both larger and smaller scales. 

Having demonstrated the broad features of the dynamo-generated field, and shown that a field 
with large-scale structure in the horizontal plane is indeed produced, it is important to seek an 
understanding of the underlying physical processes responsible for field generation. In particular, 
is it possible to distinguish between the various scenarios outlined in §[2]? To this end, in the 
following section we try to answer this question by comparing the dynamo properties of the 
actual convective flows with those of related flows obtained by the removal of selected Fourier 
modes, a process we term 'filtration'. 



5. Filtered Flows 

5.1. The Filtration Process 
The process of filtration that we employ is ess entially th a t of lo w- and high-pass filtration 



first introduced into the study of turbulence by lObukhov I (119411) : in the context of isolating 



the important modes for dynamo action, the idea of spectral filtering has been explored by 
Tobias & Cattaneo I d2008l) . 



By virtue of the periodicity in the x and y directions, all variables can be expressed as a sum 
of Fourier modes of the form 

f(z, k x , k y ,t) exp ±i{2it\~ 1 (k x x + k y y)), (5.1) 

where k x and k y are integers. If we denote a cut-off wave number by fc cut then the filtration takes 
one of the following forms: 

(a) short wavelength (SW) cutoff (i.e. long wavelengths retained): set to zero the amplitudes of 
all modes for which k 

(b) long wavelength (LW) cutoff (i.e. short wavelengths retained, plus the shear): set to zero 
the amplitudes of all modes for which k = mm(\k x \, \ k y \) < k cut , but retain the mode (0, 1) 
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Figure 10. Density plots of the vertical velocity close to the upper boundary for (a) SWC30, (c) SWC20 
and (d) SWC10. The residual small-scale flow removed by the filtration SWC30 is shown in (b). Plots (a), 
(6), and (c) are on the same scale; the residual velocity in (b) is of much smaller amplitude and is scaled 
independently. White denotes upward velocity, black falling velocity. 



corresponding to the shear. It is worth stressing that the amplitude of the (0, 1) mode emerges 
from the interaction of the convection with the imposed shear of amplitude Uq; for the flows 
considered in this section, the energy in the (0,1) mode is about 90% of the target energy. 

The filtration is applied in k x and k y since we are addressing the issues of scale separation and 
large-scale field generation in the horizontal plane. For completeness, we have also performed 
some runs in which filtration has been applied to the vertical spectrum, but we have found that 
this further filtration makes no significant difference to the properties of the magnetic fields that 
result. 

Whichever the filtration adopted, the procedure is as follows: 

(1) Solve the momentum and heat equations at full resolution; 

(2) At each time step perform the filtering to produce a filtered velocity Uf together with the 
shear; 

(3) Solve the induction equation (13.31) at full resolution with u replaced by Uf . 

It is helpful to introduce some notation: thus SWCn denotes a short-wave cutoff at k cut = n 
and LWCn a long- wave cutoff at k cut = n. In the following subsections we investigate the short- 
and long-wave cutoffs for a range of values of fc cut for the case of Ra — 150 000, Uq = 1000. 
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Figure 11. Magnetic energy versus time, with Uo = 1000 and for various shortwave cutoffs SWCn: 
(a) all modes retained, (b) n = 20, (c) n = 10, (d) n = 5,(e)n = 4, (/) n = 3, (g) n = 2. 



5.2. Short Wavelength Cutoffs 

The influence on the convection of filtration via a short wavelength cutoff is exhibited clearly in 
plots of the vertical velocity, as shown in figure [TO] As is to be expected, notable changes come 
about when k~* t is comparable with the number of convective cells across the domain. For the 
example shown in figure [lOl there is little readily appreciable difference between the SWC30 
flow and the full flow; the high frequency residual velocity that is removed by the filtration is 
shown in figure [TOtfr). 

Before discussing the influence of short wave cutoffs on flows with shear, it is important to 
understand the dynamo properties of an unsheared flow subject to the same filtration process. For 
the case of Ra = 150 000, which, recall, does not act as a dynamo in the absence of shear, the 
dynamo properties initially improve as k cut is decreased, owing to the removal of the damping 
effect of the small scales; indeed there is a range of k cut around 20 for which the filtered flow 
acts as a dynamo. However, when k cut is sufficiently small so as to exclude the energy-containing 
modes of the convection, the filtered flow becomes too feeble to support dynamo action. 

Analogous behaviour occurs in the sheared case, as shown in figure EU which plots the mag- 
netic energy versus time for various filtrations with Uo = 1000. As in the unsheared case, dy- 
namo action is enhanced slightly by the removal of just the smallest scales. In comparison with 
the unsheared case, the value of k cut at which a significant reduction in the growth rate occurs is 
now rather smaller, reflecting the importance of modes of scale intermediate between that of the 
shear and that of the original (unsheared) convection; figure [2 clearly shows these longer scales. 
The SWCn (n = 3, 4, 5) flows are still able to support dynamo action, although the growth is 
weak and somewhat irregular. It is important to note that the dynamo fails if k cut = 2, thereby 
showing that the dynamo process does not depend on the largest velocity scales alone, but must 
rely crucially on velocity scales comparable with those of the convection. 

The fact that the growth rate of, for example, the SWC20 and SWC10 dynamos is similar 
to that resulting from the full, unfiltered flow does not, of itself, mean that the same dynamo 
mode is being excited. It is, in addition, necessary to examine the morphology of the magnetic 
field generated. Comparison of figure [121 with figure [71 which show the same localization and 
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Figure 12. Density plots of B x at the upper boundary, for various shortwave cutoffs SWCn: (a) n — 20, 
(b) n — 10, (c) n — 5. Colour table as in figure [7] 
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FIGURE 13. Magnetic energy versus time for various long- wave cutoffs LWCn: (a) all modes retained, (b) 

n 20, (c) n 10, (d) n = 5,(e)n = 4, (/) n = 3,(g)n = 2. 

striated structures, does though confirm that the dynamo mechanism is identical in the filtered 
and unfiltered cases. 

5.3. Long Wavelength Cutoffs 

The idea behind implementing long wavelength cutoffs is to obtain velocity fields with a strict 
scale separation between the large scale of the shear and a much smaller convection scale; in this 
way we can find to what extent the dynamo is of a classical mean field auj type. It should first 
be noted that there is no dynamo action resulting from flows with long-wave cutoffs that also 
discard the target shear mode. Although perhaps not too surprising, this rules out the possibility 
that the influence of the shear on the small scales is such that the small scales, of themselves, 
become capable of acting as a dynamo. 

Figure [13] shows the temporal growth of the magnetic energy for LWCn flows for a range 
of values of n. For sufficiently large n, the small scales retained, in conjunction with the shear 
mode, do not support dynamo action. However, the n = 20 flow does act as a dynamo, with a 
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Figure 14. Density plots of B x at the upper boundary, for various long- wave cutoffs LWCn: (a) n = 20, 
(b) n — 10, (c) n — 5. Colour table as in figure [7] 



very well defined uniform growth of the field with time, albeit with a slow growth rate. As can be 
seen from figure[l3l decreasing n further leads to more efficient dynamo action, with the dynamo 
growth rate of the LWC5 flow becoming comparable with that of the full convective flow. It is 
surprising that dynamo action for the LWC2 flow is much weaker than that for the full flow; it is 
though possible to shed some light on this by consideration of the structure of the magnetic field 
generated. 

Simply from inspection of the growth rate, we cannot rule out the possibility that the dynamo 
mechanisms of the full flow and the LWC flows are the same; for example, the growth rates of 
LWC10 and LWC5 are very similar to that of the full flow. However in this case, examination 
of the magnetic field reveals that of the LWC dynamos to be very different in spatial structure 
from that resulting from the full flow. Figure [14] shows B x at the upper boundary for three LWC 
flows. In marked contrast to the magnetic field generated by the full flow (figure [7]), the field is 
essentially zero in A/2 < y < A and is confined to two bands in < y < A/2. It is of interest 
to note that magnetic field is expelled from regions of vigorous flow in A/2 < y < A and from 
the region of strong shear inO < y < A/2. Figure ITU exhibits a large-scale modulation in y 
in addition to small-scale fluctuations. Therefore it might be thought that such a dynamo could 
be understood within the mean field framework with an averaging that allows for y modulation. 
However, direct calculation of the a-effect reveals that a is large where the field is weak and so 
it is not clear that such a mean field description is appropriate. The anomalous behaviour of the 
LWC2 flow in figure Q3] can be understood in terms of the existence of two different types of 
dynamo mechanism. The field generated by LWC2 is of the type exhibited in figure [14] for the 
higher LWCn modes, and for which the addition of low n modes is eventually detrimental to 
dynamo action. It is only when the n = 1 modes are included that the true dynamo is recovered. 

6. Higher Rayleigh Number 

For Ra > 170 000 (when Ta = 500 000), the convective fl ow supports kinematic dyn amo 
action even in the absence of an imposed shear flow; as shown by lCattaneo & Hughes I d2006h . the 
generated field is small-scale. The standard formulation of mean field dynamo theory proceeds 
on the assumption that small-scale dynamo action is not sustainable and that any small-scale field 
results only from the interaction between a large-scale field and a small-scale velocity. Although 
here the situation is somewhat different in terms of the description of the large-scale magnetic 
field, it is nonetheless of interest to explore whether the results of § £1415 1 are critically dependent 
on the lack of dynamo action in the absence of shear, or if they are more widely applicable. Here 
we concentrate on the case of Ra = 250 000. 
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FIGURE 15. Dynamo growth rates versus Uo for Ra = 250 000 and A = 10. 

(a) (b) 





FIGURE 16. Snapshots of the temperature perturbations close to the upper boundary for Ra 
and two different values of the shear flow: (a) Uo = 600, (b) Uo = 1630. 



250 000 



Figure [T5]plots the dynamo growth rate versus Uo for Ra = 250 000. Here the higher Rayleigh 
number leads to a greater kinetic energy in the absence of the shear flow, leading to S « 
C/o/1000. The largest value of Uo shown is Uo = 1630, which has the same value of the shear 
parameter S as the flow with Uo = 1000 at Ra = 150 000. The incorporation of shear is again 
destabilising, with a similar, non- straightforward, dependence of the growth rate on Uo as exhib- 
ited at the lower Rayleigh number (cf. figure . Three different regimes can be identified. For 
Uo < 400 (S < 0.4), the shear is not sufficiently strong to change the basic convection pattern 
and hence the resulting dynamo action. For 500 < Uo < 1100, the small-scale convection is 
modulated by large-scale vorticity, as shown in figure [T6ta). For Uo > 1200, the convective 
pattern clearly reflects the influence of the target shear flow, as shown in figure [T6l b). Note that 
although figureEJc) and figure [T6Ffr) show flows with the same formal value of the shear param- 
eter S, and are indeed similar in structure, the shear at the lower Ra is slightly more dominant in 
stretching out the convective cells inO < y < A/2. 

Comparison of the growth rate dependencies and convective flow patterns for the two differ- 
ent Rayleigh numbers suggests that, once the shear flow is influential, the underlying dynamo 
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FIGURE 17. Density plots of B x at the upper boundary, for Ra — 250 000: (a) all modes retained, 
(b) SWC10, (c) LWC10. Colour table as in figure[7] 



mechanism is the same in the two cases; the fact that the flow acts as a dynamo at the higher 
Ra in the absence of shear would therefore appear immaterial. To confirm this it is though also 
necessary to look at the structure of the magnetic field generated at the higher value of Ra and 
to verify that the effects of the filtration process are similar to those discussed in §5\ As for the 
lower Rayleigh number example, we have examined the nature of the dynamo action resulting 
from short and long wavelength cutoffs of the flow. The overall trend is found to be the same. 
For the short wavelength cutoffs SWCn, the dynamo growth rate first increases as n decreases, 
is maximized at some n in the range 10 < n < 20, and then decreases rapidly for n < 5; thus a 
range of spectral modes is required for efficient dynamo action. As shown in figures [T7fq, b), the 
magnetic field generated by the short wavelength cutoffs is consistent with that generated by the 
entire flow. For long wavelength cutoffs LWCn, dynamo action ensues for n < 30 (a somewhat 
higher value of n than for the Ra = 150 000 case) and the growth rate initially increases as n 
decreases; again though, as can be seen from figure flTf c), the resulting magnetic field is of a 
very different form to that generated by the full flow. 



7. Discussion 

The research reported i n this paper has allowed u s to gain a full understanding of the phe- 
nomenon first described in Hu ghes & Proctor I d2009h . which indicated that the incorporation of 



a large-scale shear flow into rotating convection promoted the generation of large-scale magnetic 
fields. We have examined the importance for field generation of the various scales in the flow 
by considering spectrally filtered velocity fields. Our main result, which we believe to be poten- 
tially significant in terms of understanding astrophysical magnetic field generation, is that the 
observed dynamo process depends for its existence on the entire range of scales from the shear 
flow down to the scale of the convective cells. In §[2] we speculated that the introduction of shear 
might enhance the efficacy of the two-scale (mean field) dynamo process for turbulent, high Rm 
flows. Instead we see that the dynamo is produced by a completely different mechanism with 
no scale separation. It thus remains an open question as to whether the mean field ansatz is ever 
appropriate in these circumstances. 

Our first aim was to confirm the earlier result that the addition of shear to a convective flow led 
to an enhancement of dynamo action. The growth rate of the dynamo i ncreases with the magni - 
tude of the shear S. Numerical studies of forced t urbulence with shear bylYousef et al. I d2008allbb, 



together with a calculation for a simple model by lHeinemann. Mc Williams & Schekochhin I d201 



Si, 



suggest that over a range of s hear amplitudes the growth rate of such a dynamo should increase 



linearly with the shear, though lProctorl (120 12|) has shown analytically that the linear scaling can- 
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not continue to arbitrarily large values of S. In the present case, for both values of the Rayleigh 
number studied, the growth with S is far from linear. This is because the nature of the flow 
changes significantly as the shear is increased, as discussed in S I4.1I 

We considered two different values of the Rayleigh number. For Ra = 150 000, the convective 
flow is not a dynamo, and so the effect of shear is crucial. For Ra = 250 000, on the other hand, 
there is a small-scale dynamo even in the absence of shear. Adding shear to this flow enhances 
dynamo action and produces a large scale component to the magnetic field, but the morphology 
of the evolving fields is very similar when S = 0(1). This suggests the interpretation that the 
action of the shear, whether on a small-scale 'non-dynamo' or a small-scale dynamo, is very 
similar in both cases. 

Whilel Yousef et aT\ ( 2008bl) considered forced rotating turbulence with shear, the earlier paper 
( Yousef et al. Il2008a|) has no superimposed rotation, and yet the dynamo appears to function in 
a very similar manner in both cases. This led us to undertake further computations to look at the 
effects of shear on dynamo action in a non-rotating layer. However, in the absence of rotation, 
the imposed shear turns out to be readily destabilized by the convective flow and, at least with 
our target shear flow, it was not possible to attain a stable state with O(l) values of the shear 
parameter S. The same problem precluded any systematic investigation of the 'shear-current 
effect' (see, for example. iRogachevskii & Kleeorin 1120071: ISridhar & Singh 11201 0|) . 

We have tried to understand the nature of the dynamo process by considering 'filtered' flows. 
Clearly the dynamo is not much influenced by the smallest scales of flow, but removing scales 
intermediate between the shear and the turbulence has a huge effect on the form of the growing 
field, but not necessarily on the growth rate. One might expect a true mean field dynamo in the 
latter case, but in fact the largest scale of variation of the field is much smaller than the scale of 
the shear and, although a reduced model might be constructed by averaging along the direction 
of the shear flow, the resulting emf cannot be represented by a mean-field coefficient of the usual 
kind. So, paradoxically, the creation of conditions for a mean field dynamo precludes a dynamo 
of mean field type! All scales of flow except the smallest are needed to describe the dynamo 
process that is observed. To date, we have only considered a simple filtering process in Fourier 
space, isotropic in the horizontal directions. Given that anisotropy is introduced by the shear flow, 
it would be interesting to consider filtrations for which fc cut is different in the x and y directions. 
More broadly, further physical insights may be gained by employing a wavelet filtration, where 
one could filter in space as well as in scale of variation. 

There are of course other ways of combining shear and convection to produce a dynamo; 
for example, shear can be created through an Ekman layer in a rotating convecting fluid (e.g. 
Pontv. Gilbert & Sulemll200ll:lzhang. Gilbert & Zhang lEooeih : alternatively, the shear might be 



produced as a thermal wind by horizontal temperature gradients. It would be of interest to know 
what scales of motion control the appearance of the dynamo in these cases. 

The present study has investigated only the kinematic phase of the dynamo; the effects of the 
Lorentz force on the flow, which will eventually lead to equilibration, have been ignored. The 
final form of the magnetic field and the relation of the shear a mplitude to the final mag netic en- 
ergy are, however, of considerable interest. Our earlier results dHughes & Proctor ll2009h suggest 
that, somewhat surprisingly, the final field amplitude is almost unaffected by the shear provided 
that the dynamo is sufficiently vigorous. We intend to return to this question in future work. 
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